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The local density of states (LDOS) is a distribution that characterizes the effect of perturbations 
on quantum systems. Recently, it was proposed a semiclassical theory for the LDOS of chaotic 
billiards and maps. This theory predicts that the LDOS is a Breit-Wigner distribution independent 
of the perturbation strength and also gives a semiclassical expression for the LDOS witdth. Here, 
we test the validity of such an approximation in quantum maps varying the degree of chaoticity, 
the region in phase space where the perturbation is applying and the intensity of the perturbation. 
We show that for highly chaotic maps or strong perturbations the semiclassical theory of the LDOS 
is accurate to describe the quantum distribution. Moreover, the width of the LDOS is also well 
represented for its semiclassical expression in the case of mixed classical dynamics. 
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I. INTRODUCTION 

The response of quantum systems to external pertur- 
bations is a problem of paramount importance in many 
areas of physics. Many of the properties of complex 
quantum system change dramatically when the system is 
perturbed, generating fundamental phenomena as quan- 
tum phase transitions, irreversibility or dissipation. The 
present development of experimental technics in quan- 
tum complex systems make the understanding and char- 
acterization of the effect of perturbations highly desir- 
able. 

The most likely suitable magnitude to characterize the 
effects of perturbations on quantum systems is the local 
density of states (LDOS). The LDOS, also called strength 
function, was introduced by Wigner [l| to understand 
the statistical properties of the wave functions of com- 
plex quantum systems. The LDOS is the profile of an 
eigenstate of an unperturbed quantum system over the 
eigenbasis of its perturbed version. To be more specific, 
let us consider a system with a one parameter dependent 
Hamiltonian H(k) with cigcnfrccuencies ojj{k) and eigen- 
states \ipj(k)}. The LDOS of an eigenstate \ipi(k )} (that 
we call unperturbed) is given by 

Pl (^,5fc) - ^ K^ (fc)|^(fc ))| 2 5(a; - ^(fc,fc )), (1) 
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with uiij(k, fco) = LUi(ko) — ujj(k) and 5k = k — kg the per- 
turbation strength. Eq. Q]shows that the LDOS is a den- 
sity of states in which the delta functions are weighed by 
the overlaps between perturbed and unperturbed states. 
In addition, the LDOS width gives an estimation of how 
many perturbed states contribute to an unperturbed one. 
Besides, it is the Fourier transform of the fidelity ampli- 
tude (FA) of the state \ip%{ko)), 

Pi (cj,5k) = T[(Mh)\e lHik)t/h e- lH ^ h \Mko))]. (2) 
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Both the FA and its absolute square value, called the 
Loschmidt echo, are important measures of sensitivity to 
perturbations and irreversibility of quantum evolutions 

The LDOS has been considered in many contexts. In 
a seminal paper, Wigner studied the LDOS in a sim- 
ple model of banded random matrices [l[ . Subsequently, 
many authors have used the LDOS to characterize the 
structure of the eigenstates of different random matrix 
models 0-0 • The LDOS has also been studied in several 
microscopic systems as for example in a Ce atom (Toj . in 
chaotic billiards [ll| or a system of a particle that evolves 
in a smooth Hamiltonian [12] ■ In addition, the LDOS has 
been studied to characterize the effect of perturbations 
in the operation of quantum computers in the presence of 
static imperfections [HI, [l4[ . It was shown that depend- 
ing on the characteristics of the system, the LDOS has 
many regimes as a function of the perturbation strength 
5k. However, all the mentioned studies have revealed a 
region of perturbation strength in which the LDOS has 
a lorentzian shape, that has been usually called Breit- 
Wigner distribution. 

A step forward has been recently made in the under- 
standing of the LDOS for chaotic systems Its rela- 
tion with the FA has been exploited to develop a semiclas- 
sical theory of the LDOS for locally perturbed billiards 
or maps, that is, when the perturbation is concentrated 
in a small region of the phase space accessible for the sys- 
tem. It was shown that the LDOS has a Lorentzian shape 
under very general perturbations of arbitrarily high in- 
tensity an a semiclassical expression for its width was 
derived. This expression only depends on the perturba- 
tion, while the properties of the system are taken into 
account through a uniform measure in phase space. The 
same results were obtained in a subsequent publication 
for maps that are globally perturbed but the dynamics 
was assumed to be completely random (l6| . 

The aim of our study is to test the validity of the semi- 
classical theory of Ref. [Tol [l6j in quantum maps when 
the perturbation is applied in all the phase space and the 
dynamics of the classical map is not completely random. 
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We also consider perturbations that act in an region of 
the phase space. We study the behavior of the LDOS 
for maps with different degree of chaoticity and intensity 
of the perturbation. For this purpose we consider two of 
the most paradigmatic systems of quantum chaos studies: 
the perturbed cat map and the Harper map. We show 
that the semiclassical approximation of the width of the 
LDOS works very well even for systems with mixed dy- 
namics in which chaos coexist with regular islands. The 
prediction of Lorenzian shape of the LDOS is fulfilled for 
highly chaotic maps or when the intensity of the pertur- 
bation is big enough. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce the dynamical systems that we have used for the 
numerical study, the cat and the Harper maps. In this 
section, we describe the main characteristics of the clas- 
sical and quantum dynamics of the maps. Sec. Mil is 
devoted to present the semiclassical theory of the LDOS 
[lH The starting point of this theory is a semi- 

classical approximation of the fidelity amplitude called 
Dephasing Representation [lj|. In Sec. IIVI we study the 
behavior of the LDOS for the systems in several situ- 
ations and test the validity of the semiclassical theory. 
We consider various degree of chaoticity and intensities 
of the perturbation. We also compare the cases ol local 
and global perturbations. Finally, we conclude with a 
summary of our results and some final remarks in Sec. 

El 



II. SYSTEMS: MAPS ON A TORUS 

An usual procedure to understand a complex behav- 
ior is to consider very simple systems in which such a 
phenomena is observed. The most simple dynamical sys- 
tems which develops all types of complexity are abstract 
maps. Due to their simplicity, classical and quantum 
maps have been very important in the development of 
classical and quantum chaos j20l - |22l |. Furthermore, many 
quantum maps have been implemented experimentally in 
previous studies p3l - l25j . 

In this paper we have used maps acting on a torus 
phase space of area A. = 1. In particular we have con- 
sidered the well known cat and Harper maps. These 
maps possess all the essential ingredients of chaotic and 
mixed dynamics and are extremely simple from a numer- 
ical point of view. 

The cat maps are linear automorphisms of the torus 
that exhibit hard chaos. Anosov's theorem [26| estab- 
lishes that the cat maps are structurally stable, that is, 
the orbits of a slightly perturbed map are conjugated to 
those of the unperturbed map by a homcomorphism. A 
perturbation of a cat map can be represented by matrices 
acting on the coordinates 



= G 



where G is a 2 x 2 matrix with integer elements chooscn 
that Tr (G) > 2 and det(G) = 1 since the maps are 
hyperbolic and conservative. We consider a perturbation 



e(q,k) = (fc/27r)[cos(27rq) - cos(47rg)], 



(4) 



with the perturbation strength k < 0.11 to satisfy the 
Anosov theorem [26|, [27| ■ To take into account different 
degrees of chaoticity, in this paper we have considered 
the following matrices G, 



G l = 



2 1 
1 1 



Go = 



80 1 
6399 80 



The corresponding Lyapunov exponents, which deter- 
mine the rate of exponential divergence of classical tra- 
jectories are Ai ~ 0.96 and A2 ~ 5.07. We note that A is 
approximately uniform over the whole phase space and 
nearly independent of k t 28| . 

Perturbed cat maps do not capture all the possible mo- 
tions of Hamiltonian systems. The most common situa- 
tion is a mixture of regular islands interspersed by chaotic 
regions. To consider this general situation the model that 
we have chosen to study is the Harper map in the unit 
square (2^ |. 



q = q — k sin 2irp 
p = p + k sin 2irq 



(mod 1), 
(mod 1), 



(5) 



where k is a parameter that controls the behavior of the 
system. This map can be understood as the stroboscopic 
version of the flow corresponding to the (kicked) Hamil- 
tonian 

H(p, q,t) = -— cos(2ttp) - — cos(27rg) ^ S(t - nk). 



e(q,k) (modi), (3) 



(6) 

This is an approximated Hamiltonian for the motion of 
an an electron in a crystal under the action of an external 
field. 

The Harper map presents a mixed dynamics that de- 
pends on the parameter k. Fig.[T]shows some phase space 
pictures for this model as an example of the underlying 
classical dynamics. As can be seen in Fig. Q] (left panel), 
the system presents a mixed dynamics with regions of 
regularity around the origin and the corners coexisting 
with chaos, in agreement with the KAM theorem [2r| . 
When the parameter k = 200 the size of the island are 
so small that it is not possible to observe without a finer 
resolution [see Fig. [T] (right panel)] . 

The quantization on the torus implies that the wave 
function should be periodic in both position and mo- 
mentum representation. If in the coordinate and mo- 
mentum representation the wave function has a period 1 
with spacing 1/N, it follows that 1 = 2irhN . Then, we 
have a Hilbert space of N dimension for a fixed value of K. 
As N takes increasing values, we reach the semiclassical 
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FIG. 1. Classical phase space of the Harper map for k — 0.3 
(left panel) and k = 200 (right panel). See text for details. 



The inverse Fourier transform of Eq. [9] the so called 
average fidelity amplitude (AFA), is the starting point of 
the a semiclassical approximation of the LDOS, 



(k)t/h-iH (k )t/h 



To evaluate Eq. [10] we have used the so called dephas- 
ing representation, a semiclassical formulation for fidelity 
amplitude which avoids the usual traj ectory-search prob- 
lem of the standard semiclassics [Hj. One of the forms 
of the FA obtained using the dephasing representation is 



limit. The position basis {ft}^ 1 (with qi — i/N) and 
momentum basis {pt} 1 ^ 1 (with pi — i/N) are related by 
the discrete Fourier transform. In this setting a quantum 
map is simply a unitary U acting on an N dimensional 
Hilbert space and evolution after n steps is given by U n . 

There is no general method for map quantization. For 
the perturbed cat map we have considered the quantiza- 
tion based on the classical propagator of Rcf. [20, H3] • In 
this case, the matrix elements of the propagator in the 
position basis are 




{gnq 2 - 2q'q + .g 22 g' : 



\Sin(27rg) - - sin(47rg)) 



(7) 



where g- t j arc the elements of the matrix G and we have 
used g\2 = 1. 

For the Harper map [29l | , the matrix elements of evolu- 
tion operator in the mixed basis of position and momenta 
are 



Ufr (q p) = e iNkcos ( 2 ^l) e iNkcos(2Trp) 



(8) 



III. SEMICLASSICAL THEORY OF THE LDOS 
OF CHAOTIC MAPS 

The LDOS as defined in Eq. Q] depends on the charac- 
teristics of the state \ipi{k )). To avoid any dependence 
with some particular characteristics of this state an av- 
erage over unperturbed states is performed. Due to the 
finite number of states in quantum maps, we average over 
all the Hilbert space. Thus, the averaged LDOS p(ui,dk) 



N 



p(uj,Sk) = — ^2pi(u},5k). 



(9) 



Of(t,6k)= W 4> (q,p)e~ lAS ^ Sk ^ h dqdp, (11) 



where ASt{q,p, Sk) is the action difference evaluated 
along the umperturbed orbit starting at (q,p) that 
evolves at a time t and W<p(q,p) is the Wigner function 
of the initial state Id). Then, 



0(t,Sk)= W(q,p)e- lASt ^ p ' 5k)/h dqdp, (12) 



where W{q,p) = {l/N)Y,Wi(q,p), with Wi(q,p) being 
the Wigner function of \ipi(ko)). For chaotic systems, 
the mean value of the Wigner function for a base of 
cigenstates is approximately a uniform distribution so 
W(q,p) = 1/V where V is the volume of the phase space. 
Therefore, 



0(t,Sk) 



1 

V 



-iAS t {q,p,Sk)/h 



dqdp. 



(13) 



Time is discrete in maps, so from now on we use the 
integer n to count time steps and V is the area of the 
phase space that in our case is equal to unity. 

In order to solve Eq. [13] for maps we need to assume 
that trajectories become uncorrelated between two suc- 
cessive hits in the perturbed region. This approximation 
is valid when the perturbation acts on an infinitesimal 
portion of phase space [l^, [13, Ell or ^ the unperturbed 
dynamics of the system is completely random [16l | . 

Here we have considered the second case, the A — > oo 
limit, by assuming that the dynamics is purely random. 
This evolution is completely stochastic in the sense that 
there is no correlation for the differ ent times of the evo- 
lution. Then, to compute 0(n,dk), we have divided the 
phase space in N c cells. The probability to jump from 
cell to any other in phase space is uniform. Therefore it 
is straightforward to show that the mean FA results 



0(n,5k) =£-E e H(AS, 1+ ...+A5*,)/A] 

jl in 



E' 



(-iASj/h) 



(14) 
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where ASjp is the action difference evaluated in the cell 
j at time p. The continuous limit is approached when 
N c —> oo resulting in 



0(n,Sk) 



,-iA.S(q,p,8k)/h 



dqdp 



(15) 



where AS(q,p, Sk) is the action difference after one step 
of the map. 

The exponential decay of Ec. [TSlcan be rewritten as 



0(n, Sk) = e 



-Tn-\-iipn 



(16) 



with 



and 



ln(| J e- iAS ^' p ' sk ^ h dqdp\). (17) 



to 
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FIG. 2. (Color online) . Witdh cr vs. T for a periodized Lorentzian 
function [Eq. I20I . The limit of a for r — > oo which corresponds for 
a constant LDOS is also plotted with red dotted line. 



tp = arg( J e- lAS ( q ' p ' 5k)/h dqdp). 



(18) 



We note that T and tp depend on the perturbation 
strength Sk. 

Now, we obtain the semiclassical expression for the 
average LDOS by the inverse Fourier transform of Eq. 

us 



distribution. In our case we are going to take the dis- 
tance about the average value of the LDOS that contains 
70 % of the probability. That is, 



where 



p{ui, Sk)diu = 0.7 



(21) 



p sc {uj,Sk) = T^(u),T,y) 



n[{uj-tp) 2 + T 2 ] 



(19) 



The phase tp determines the location of the center of the 
Lorenzian function and T its width. 

Finally, we have to take into account the fact that the 
spectrum of a map is periodic because of a compact phase 
space. This periodicity changes the form of the LDOS 
into a periodized Lorentzian function 



= E 



7r[(w - tp - 2nj) 2 + T 2 ] ' 



(20) 



The same semiclassical expressions for the LDOS were 
obtained in Ref. [HI, G3 when the perturbation acts in 
a region of the phase space of area a 0. 

A magnitude that has physical interest is the width a 
of the LDOS which is a measure of the number of per- 
turbed states that are needed to describe a unperturbed 
one. Therefore, this quantity offers clear information 
about the effect of perturbations on on a quantum sys- 
tem. Moreover, the width of the LDOS determines for 
some regime of the perturbation, the rate of fidelity decay 
under imperfect motion reversal (the Loschmidt echo). 
There are different ways of determining this width of a 



(u) = / wp(ui,Sk)du> 



(22) 



We show in Fig. [2] the relation between its width a 
and r for the periodized Lorentzian function of Eq. [201 ■ 



IV. RESULTS 

The main interest of a semiclassical theory is to de- 
scribe quantum mechanical quantities using classical in- 
formation. In this section we show the behavior of the 
LDOS for the quantum maps presented before and test 
the validity of the semiclassical approximation of the 
LDOS described in the previous section. The aim this 
section is to compare the approximated p sc and a sc with 
the corresponding exact quantum values. The latter are 
numerically computed by diagonalization of the evolution 
operators of Eq. [7] and [H 

The semiclassical approximation of the LDOS is com- 
pletely determined by T and tp that are obtained with the 
calculation of the integral of Eq. [15] To avoid the de- 
pendence of the results with the dimension of the Hilbert 
space N we have considered all the studied quantities as 
a function of the scaled strength of the perturbation 



X=(k- k )/{2irh) = SkN. 



(23) 



5 



71 


9 

r 

6 







i ■ 




1 1 


1 1 1 1 1 1 






i i i i i i i i i i 


— 1 

'-) 


u 


b 

i 


u 


1 1 

\) 

1 


1,1,1, 



30 60 90 120 150 



X 




X 

FIG. 4. (Color online) Width a of the LDOS as a function of 
the scaled perturbation strength X = — ko)N for the cat map 
with Gi (O) and G2 (X). The red solid line is the semiclassical 
approximation of a. The number of states of the Hilbert space iV = 
2000 and fco = 0.01. We indicate with arrows the perturbations 
strength of the LDOS displayed in Fig. 



FIG. 3. (Color online) Y and the phase ip as a function of the 
scaled pertubation \ for the perturbed cat map. 



In all the calculations included in this section the number 
of states of the Hilbert space set as N = 2000. 



A. Peturbed cat map 

The action difference for one iteration of the perturbed 
cat map described in Sec. [TT| is given by 



AS(q,p,6k) 



5k 
4^2 



sin(27rg) — — sin(47rg) 



(24) 



Using Eg. [24l.[T7land[T8lwe compute V and p>. In Fig. [3] 
we plot r and <p for the perturbed cat map as a function 
of the scaled perturbation strength x- We can see that 
for perturbation of Eq. [4] ip has only two possible values 
cither or 71". 

We firstly compare the semiclassical approximation of 
the width of the LDOS with the corresponding quantum 
value. For this reason the width of the LDOS has been 
computed for the cat map using k = 0.01 to avoid all 
the arithmetic peculiarities of the cat map (k = 0), which 
account for the non-generic spectral statistics [30j. In 
Fig. [4] the width of the LDOS is shown for the cat maps 
with Gi and Gi- The semiclassical approximation er sc , 
plotted in solid line, works extremely well for both cat 
maps in the whole range of considered perturbations. 

The width of the LDOS a for the cat maps has two 
clearly different regimes [ Fig. [3] . For small perturba- 
tion strength when % < 10 it presents a quadratic be- 
havior that is usually called Fermi Golden Rule regime. 



Conversely for greater strength the width is an oscillat- 
ing function. In order to understand the behavior of a sc 
when % — >■ 00 we have used the stationary phase approx- 
imation method to solve the integral Eq. [T3] obtaining 



r 



log[l/Vx] for x 



therefore the width a sc — > 0.7ir that corresponds to a 
uniform distribution. 

At this point we would like to see how good the 
semiclassical approximation of the LDOS can describe 
the complete distribution. We have therefore computdc 
the quantum LDOS for several values of perturbation 
strength for the cat maps Gi and G2. In Fig. [5] we com- 
pare the LDOS with its semiclassical approximation for 
perturbations indicated in Fig. 2]with arrows. Fig. [5{a) 
corresponds to the most chaotic case Gi- In the main 
plot x = 159.6, in the left inset x = 14-4 and in the right 
inset x — 18. We can see that the semiclassical approx- 
imation works very well for all the perturbations, that 
is, the LDOS is a periodized Lorentzian function indis- 
tinctly of the perturbation strength. Left and right inset 
corresponds to approximately the same width of the dis- 
tribution but in the left figure (p = and ip = n for the 
right so in this case the periodized Lorentzian is centered 
in lu = it. As can be seen in Fig. [3j near x ~ 15, the 
phase <p has a discontinuity and jumps from to ir, for 
this reason the center of the LDOS changes from to = 
to uj = ir. Similar behavior occurs in the other disconti- 
nuities of ip near x ~ 50 and 70. 

In Fig. [5jb) the results for the cat map with the ma- 
trix Gi are shown. In the figure main panel, we show 
that for big perturbation strength after the quadratic 
regime (x = 159.6) the LDOS is well described by 
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FIG. 5. (Color online) LDOS p (points) and its semiclassical ap- 
proximation (solid line), (a) cat map of Eq. [3] with G = G2 for a 
scaled perturbation strengh x = 159.6 (main plot), \ = 14.4 (left 
inset) and \ = 18.0 (right inset), (b) cat map of Eq. [3]with G = G\ 
for x = 159.6, x = 14-4 (left inset) and x = 18-0 (right inset). 



the semiclassical Lorentzian distribution. Conversely for 
smaller perturbations strength the LDOS does not show 
a Lorentzian behavior [see inset of Fig. Efb)] . To under- 
stand this behavior we show in Fig. [5] the mean value 
of the fidelity amplitude 0(n) for \ = 14-4 of both cat 
maps with G\ (o) and G% (x) which corresponds to the 
inverse Fourier transform of the LDOS plotted in the left 
inset of Fig. [5{a) and (b). We see that in the case in 
which the LDOS is not a periodized Lorentzian function 
the corresponding 0(n) has a big revival (at n = 4). This 
kind of behavior, known as survival collapse after which 
the largest revivals appear, was observed in a spin chain 
[3^ and can be the cause for non-Markovian quantum 
evolutions (33|. 

We test now the validity of the semiclassical approxi- 
mation of the LDOS for local perturbation. For this rea- 
son the perturbation is applied in a q strip from qo = 0.25 
to qi = 0.46 so the area of the perturbed region is 
a = AqAp = qi — qo = 0.21. In Fig. [7| we show T 
and ip as a function of the scaled perturbation strength 
£ computed using Eq. [17] and [TBI We can see that for 
this local perturbation ip is an oscillating function so the 
semiclassical approximation of the LDOS is periodized 
Lorentzian function with an oscillating mean value. In 
Fig- Eltop panel) the mean value of the exact LDOS it 
is also plotted with (□) showing that the semiclassical ip 
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FIG. 6. (Color online) Mean value of the amplitude fidelity 
0(n) for the cat map with Gi (o) and G2 (X) with pertur- 
bation strength \ = 14.4. The exponential decay given by 
exp(— Tn) is plotted with solid blue line. We also plot with 
dotted red line the exponential decay given by exp( — An/2) 
with A the Lyapunov exponent for the cat map with Gi. 




30 60 90 120 150 



% 

FIG. 7. (Color online) T and the phase 93 as a function of the scaled 
pertubation x for the cat map when the perturbation is applied in 
a q strip from qo = 0.25 to qi = 0.46. The mean value of the 
quantum LDOS is also plotted with (□). 



describe very well this quantity. 

The LDOS is also very well approximated by the semi- 
classical LDOS for all the perturbations strength that we 
have studied. In Fig. [8jwc show the LDOS for in \ = 8 
when the width grows cuadratically (FGR regime) and 
for x = 28 when the width shows an oscillating behav- 
ior. The semiclassical approximation is plotted with solid 
line. In the inset of Fig[5]wc show the width of the LDOS 
for this local perturbation and its semiclassical approxi- 
mation. We can clearly see that the a sc works very well 
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FIG. 8. (Color online) LDOS p for a local perturbations. The 
cat map with Gi is perturbed from go = 0.25 to q\ = 0.46. The 
scaled perturbation strength \ = 8 (A) and \ = 28 (□). The 
semiclassical approximation of the LDOS is plotted with red solid 
line. Inset: Width a of the LDOS as a function of the scaled 
perturbation strength \ (°) an d with red solid line the semiclassical 
approximation. 



for local perturbation. It is noteworthy that all the calcu- 
lations for local perturbations were done using the map 
G\ showing that when the perturbation is applying is a 
small region of the phase space less degree of chaoticity 
is needed for the semiclassical LDOS to be accurate. 



B. Harper map 

We have studied the LDOS of the Harper map using 
the evolution operator of Eq. [8] with k = k + 8k. The pa- 
rameter Sk is the perturbation strength and as we have 
used for the cat map, the scaled perturbation straight 
X = SkN. In this case the action difference for one iter- 
ation of the Harper map is given by 



AS(q,p,Sk) 



5k\ 

— [cos(2ttp) + cos(2tt</)] . (25) 
2n J 



where q' is given by Eq. 

We have considered as unperturbed system the cases 
with fco = 0.30 [mixed dynamics, Fig. [JJleft panel)] and 
ko = 200 [chaotic dynamics, Fig. QJright panel)]. Using 
Eq. ITTl [P8l andl20lwc compute T, ip and the corresponding 
semiclassical approximation of the LDOS. In Fig. [5] we 
show r as a function of the scaled perturbation strength 
X- For the action difference of the Harper map [Eq. [25] 
we have obtained that ip = 0. 

In Fig. [10] we show the width of the LDOS for the 
Harper map and the corresponding semiclassical approx- 
imation. When the dynamics of the Harper map is com- 
pletely chaotic, the semiclassical a sc works well as ex- 
pected. Surprisingly, the semiclassical approximation 
works reasonably well even for mixed dynamics. This 




FIG. 9. (Color online) T as a function of the scaled pertubation 
strength \ f° r the Harper map. 




FIG. 10. (Color online) Width a of the LDOS as a function of 
the scaled perturbation strength \ f° r the Harper map with 
fco = 0.3 (o) and fc = 200 (X). In red line it is plotted the 
semiclassical approximation a 3C . We indicate with arrows the 
perturbations strength of the LDOS displayed in Fig. [11] 



agreement is more noticeable for bigger x- The explana- 
tion of this unexpected behavior is as follows. Eq. [TBI is 
exact for one time step (n = 1) and if the perturbation 
strength is big enough the fidelity amplitude decays in 
this short time. Therefore, this short time decay gives 
the width of the Fourier transform which is the LDOS. 

In Fig. [TTJ we show the LDOS for the Harper map. 
Although the semiclassical width of the LDOS a sc works 
well for mixed dynamics, the complete distribution is not 
well reproduced by a periodized Lorentzian distribution. 
This is shown in the inset of Fig. [TTJ (b) for the Harper 
map with fco = 0.3 and x = If the perturbation 

strength is bigger [Fig. [TTJ (a) (main plot)] the semi- 
classical theory works reasonably well but the quantum 
LDOS is a more fluctuating function than the chaotic 
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FIG. 11. (Color online) p (points) and its semiclassical ap- 
proximation p ac (solid line) for the Harper map. (a) fco = 200. 
In the main plot x — 17 an d in the inset \ = 1-8. In both 
plot it is seen that the semiclassical LDOS describes the full 
quantum result, (b) fco = 0.3 (mixed dynamics). In the main 
plot x = 17 and in the inset x — 1-8. 



case [see Fig. [IT] (a)] . As expected, the semiclassical 
theory works well for the case of fco = 200 in which the 
Harper map is fully chaotic (Fig. ITTI (a)). 



V. CONCLUSIONS 

The reaction of a system to perturbations is a fun- 
damental problem in quantum mechanics. In this pa- 
per we have made a detailed analysis of the response to 
perturbations of the simplest quantum systems, which 
can have complex classical dynamics. For this reason, 
we have studied the LDOS in the perturbed cat map, a 
completely chaotic system and the Harper map which has 
mixed dynamics. Our fundamental goal was to discuss 
the validity of a semiclassical theory of LDOS that has 



been recently developed [HI, [l6| . This theory is based on 
the relation of the LDOS with the fidelity amplitude, a 
measure of irreversibility and sensitivity to perturbations 
of quantum systems. Furthermore, it uses the dephasing 
representation of the fidelity amplitude, a semiclassical 
formulation that avoids the usual problems of semiclas- 
sical theories. The main assumption of the semiclassical 
theory of LDOS is that the trajectories get uncorrclatcd 
after one step of the map. This condition is fullfilcd if the 
dynamics is completely random or when the perturbation 
is applied in a infinitesimal region of the phase space. 
Due to the fact that these conditions are not achieved in 
dynamical system, we tested the validity of such a semi- 
classical theory of the LDOS. 

We have analyzed various situations: local and global 
perturbations and also we have varied the degree of 
chaoticity. We show that the LDOS is very well described 
by its semiclassical expression when the map is highly 
chaotic, either if the perturbation is localized in phase 
space, or when the perturbation strength is big enough 
. We remark that in these cases the semiclassical LDOS 
completely reproduces the quantum version without any 
fiting parameters. Wc have studied the case of mixed 
dynamics and surprisingly enough our results show that 
the semiclassical width of the LDOS describes the full 
quantum version even in this case. 

We would like to highlite that our results could be of 
importance in the study of the LDOS of billiards. In- 
deed, the behavior of a billiard system has many resem- 
blances with maps. For example, the classical dynamics 
of a billiard can be described by a map on the bound- 
ary. Quantum billiards are realistic systems that can be 
constructed in experimental setups of several nature. In 
fact, there are cavities of microwave, acoustic or optical 
wave. The semiclassical approximation of the width of 
the LDOS has been successfully applied in billiard that 
has been perturbed both locally [H|,[3l| and globally [l5j |. 
However, in these works the behavior of the whole dis- 
tribution was not properly discussed. Further insight on 
the LDOS of this systems will be part of future studies. 
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